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I Abstract: Some connections between operator theory and wavelet anal- 

fH , ysis: Since the mid eighties, it has become clear that key tools in wavelet 

analysis rely crucially on operator theory. While isolated variations of 

a wavelets, and wavelet constructions had previously been known, since 
Haar in 1910, it was the advent of multiresolutions, and sub-band fil- 
tering techniques which provided the tools for our ability to now eas- 
ily create efficient algorithms, ready for a rich variety of applications to 
' practical tasks. Part of the underpinning for this development in wavelet 

\—{ ' analysis is operator theory. This will be presented in the lectures, and 

we will also point to a number of developments in operator theory which 
in turn derive from wavelet problems, but which are of independent in- 
terest in mathematics. Some of the material will build on chapters in 
a new wavelet book, co-authored by the speaker and Ola Bratteli, see 
http:/ /www. math. uiowa. edu/ - jorgen/ . 
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One cannot expect any serious understanding of what wavelet analysis 
means without a deep knowledge of the corresponding operator theory. 

— Yves Meyer* 

1. Introduction 

While this series of four lectures will be on the subject of wavelets, the 
emphasis will be on some interconnections between topics in the mathe- 
matics of wavelets and other areas, both within mathematics and outside. 
Connections to operator theory, to quantum theory, and especially to sig- 
nal processing will be studied. Concepts such as high-pass and low-pass 
filters have become synonymous with wavelet tools, but they have also had 
a significance from the very start of signal processing, for example early 
telephone signals over transatlantic cables. This was long before the much 
more recent advances in wavelets which started in the mid-1980's (as a re- 
sumption, in fact, of ideas going back to Alfred Haar [HaalO] much earlier). 
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1.1. Index of terminology in math and in engineering 

Since the inid-1980's wavelet mathematics has served to some extent as a 
clearing house for ideas from diverse areas from mathematics, from engi- 
neering, as well as from other areas of science, such as quantum theory 
and optics. This makes the interdisciplinary communication difficult, as 
the lingo differs from field to field; even to the degree that the same term 
might have a different name to some wavelet practitioners from what is 
has to others. In recognition of this fact, Chapter 1 in the recent wavelet 
book [Br Jo02b] samples a little dictionary of relevant terms. Parts of it are 
reproduced here: 

Terminology 

• multiresolution: — real world: a set of band-pass-filtered com- 
ponent images, assembled into a mosaic of resolution bands, each 
resolution tied to a finer one and a coarser one. 

— mathematics: used in wavelet analysis and fractal analysis, mul- 
tiresolutions arc systems of closed subspaces in a Hilbert space, 
such as (M), with the subspaces nested, each subspace repre- 
senting a resolution, and the relative complement subspaces rep- 
resenting the detail which is added in getting to the next finer 
resolution subspace. 

• matrix function: a function from the circle, or the one-torus, 
taking values in a group of N-hy-N complex matrices. 

• wavelet: a function i/j, or a finite system of fmictions {ijji}. such 
that for some scale number N and a lattice of translation points on 
M, say Z, a basis for (M) can be built consisting of the functions 
NitlJi{N^x-k), j,kGZ. 

Then dulcet m,usic swelled 

Concordant with the life-strings of the soul; 

It throbbed in sweet and languid beatings there, 

Catching new life from transitory death; 

Like the vague sighings of a wind at even 

That wakes the wavelets of the slumbering sea. . . 

— Shelley, Queen Mab 

• subband filter: — engineering: signals are viewed as functions of 
time and frequency, the frequency function resulting from a trans- 
form of the time function; the frequency variable is broken up into 
bands, and up-sampling and down-sampling are combined with a 



February 1, 2008 



Master Review Vol. 9in x 6in - (for Lecture Note Series, IMS, NUS) 



umfwaspw 



Unitary Matrix Functions, Algorithms, Wavelets 



5 



filtering of tlie frequencies in making the connection from one band 
to the next. 

— wavelets: scaUng is used in passing from one resolution V to 
the next; if a scale N is used from V to the next finer resolution, 
then scaling by jf takes F to a coarser resolution Vi represented 
by a subspace of V, but there is a set of functions which serve as 
multipliers when relating V to Vi, and they are called subband 
filters. 

• cascades: — real world: a system of successive refinements which 
pass from a scale to a finer one, and so on; used for example 

in graphics algorithms: starting with control points, a refinement 
matrix and masking coefficients are used in a cascade algorithm 
yielding a cascade of masking points and a cascade approximation 
to a picture. 

— wavelets: in one dimension the scaling is by a number and a 



the initial step for the cascades; when the masking coefficients are 
chosen the cascade approximation leads to a scaling function. 
• scaling function: a function, or a distribution, (p, defined on the 
real line R which has the property that, for some integer A'^ > 1, 



of the linear span of the set of translated functions . . . ,ip{x + 1), 
ip{x), ip{x-l),ip{x-2),.... 

• logic gates: — in computation the classical logic gates are real- 
ized as computers, for example as electronic switching circuits with 
two-level voltages, say high and low. Several gates have two input 
voltages and one output, each one allowing switching between high 
and low: The output of the AND gate is high if and only if both 
inputs are high. The XOR gate has high output if and only if one 
of the inputs, but not more than one, is high. 

• qubits: — in physics and in computation: qubits arc the quantum 
analogue of the classical bits and 1 which are the letters of classi- 
cal computers, the qubits are formed of two-level quantum systems, 
electrons in a magnetic field or polarized photons, and they are rep- 
resented in Dirac's formalism |0) and |1); quantum theory allows 
superpositions, so states = a |0) 10), a, 6 G C, |a|^ -I- \b\^ = 1, 
are also admitted, and computation in the quantum realm allows 
a continuum of states, as opposed to just the two classical bits. 
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— mathematics: a chosen and distinguished basis for the two- 
dimensional Hilbert space consisting of orthogonal unit vectors, 
denoted |0), 

• universality: — classical computing: the property of a set of logic 
gates that they suffice for the implementation of every program; or 
of a single gate that, taken together with the NOT gate, it suffices 
for the implementation of every program. 

— quantum computing: the property of a set S of basic quantum 
gates that every (invertible) gate can be written as a sequence of 
steps using only gates from S. Usually S may be chosen to consist 
of one-qubit gates and a distinguished tensor gate t. An example of 
a choice for t is CNOT. An alternative universal one is the Toffoli 
gate. 

— m,athematics: the property of a set S of basic unitary matrices 
that for every n and every u G U2" (C), there is a factorization 
u = siS2---Sk, Si G S, with the understanding that the factors 
Si are inserted in a chosen tensor configuration of the quantum 
register • • • C^. Note that the factors Sj, the number k, and 

n times 

the configuration of the Sj's all depend on n and the gate u e 
U2" (C) to be studied. The quantum wavelet algorithm (2.2.6) is 
an example of such a matrix u. 

• chaos: a small variation or disturbance in the initial states or input 
of some system giving rise to a disproportionate, or exponentially 
growing, deviation in the resulting output trajectory, or output 
data. The term is used more generally, denoting rather drastic 
forms of instability; and it is measured by the use of statistical 
devices, or averaging methods. 

• GLat (C): the general linear group of all complex N x N invertible 
matrices. 

• Ujv (C): = { A e GLiv (C) I AA* = Icn } where A* denotes the ad- 
joint matrix, i.e., (^4*)^^ = Aj,i. 

• transfer operator (transition operator): — in probability: An 
operator which transforms signals s from input Sin to output Sout • 
The signals are represented as functions on some set E. In the sim- 
plest case, the operator is linear and given in terms of conditional 
probabilities p {x, y). The number p {x, y) may represent the prob- 
ability of a transition from y to x where x and y are points in the 
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set E. Then 

Sout {x) = ^p {x, y) Sin iv) ■ 

— in computation: Let X and Y be functions on a set E, both 
taking values in {0, 1}. Let Y be the initial state of the bit, and 
X the final state of the bit. If the process is governed by a prob- 
ability distribution P, then the transition probabilities p{x,y) := 
P{{X = x\ Y = y}) are conditional probabilities: i.e., p {x, y) is 
the probability of a final bit value x given an initial value y, and 
we have 

p {{X = x}) = Y,p y) P {{y = y}) ■ 

— in wavelet theory: Let N G Z_|_, and let be a positive function 
on T = { e C I = 1 }, for example W = |mo|^ where mo is 
some low-pass wavelet filter with N bands. (Positivity is only in 
the sense W >0, nonnegative, and the function W may vanish on 
a subset of T.) Then define a function p on T x T as follows: 

\ for all other values of w. 

We arrive at the transfer operator Rw, i.e., the operator trans- 
forming functions on T as follows: 

Sout (z) = (RwSin) i^) = ^ X! ^ ('^) *™ • 

w'^ =z 

• coherence: — in mathematics and physics: The vectors ipi that 
make up a tight frame, one which is not an orthonormal basis, 
are said to be subjected to coherence. So coherent vector systems 
in Hilbert space are viewed as bases which generalize the more 
standard concept of orthonormal bases from harmonic analysis. A 
striking feature of the wavelets with compact support, which are 
based on scaling, is that the varieties of the two kinds of bases can 
be well understood geometrically. For example, the collapse of the 
wavelet orthogonality relations, degenerating into coherent vectors, 
happens on a subvariety of a lower dimension. 
More generally, coherent vectors in mathematical physics often 
arise with a continuous index, even if the Hilbert space is sepa- 
rable, i.e., has a countable orthonormal basis. This is illustrated by 
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a vector system {tl>r,s}j which should be thought of as a continuous 
analogue, i.e., a version where a sum gets replaced with an integral 



For more details, see also Section 3.3 of [Dau92] and Chapter 3 of 
[Kai94]. 

In quantum mechanics, one talks, for example, about coherent 
states in connection with wavefunctions of the harmonic oscilla- 
tor. Combinations of stationary wavefunctions from different en- 
ergy eigenvalues vary periodically in time, and the question is which 
of the continuously varying wavefunctions one may use to expand 
an unknown function in without encountering ovcrcompleteness of 
the basis. The methods of "coherent states" are methods for us- 
ing these kinds of functions (which fit some problems elegantly) 
while avoiding the difficulties of ovcrcompleteness. The term "co- 
herent" applies when you succeed in avoiding those difficulties by 
some means or other. Of course, for students who have just learned 
about the classic complete orthonormal basis of stationary eigen- 
functions, "coherent state" methods at first may seem like a daring 
relaxation of the rules of orthogonality, so that the term seems to 
stand for total freedom! 

1.1.1. Some background on Hilbert space 

Wavelet theory is the art of finding a special kind of basis in Hilbert space. 

Let H he a Hilbert space over C and denote the inner product ( • | • ). For 
us, it is assumed linear in the second variable. liTL = L"^ (R), then 





(1.1.1) 



lin = e (Z), then 




(1.1.2) 



Let T = R/27rZ. If H = 



(T), then 




(1.1.3) 
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Functions f & L"^ (T) have Fourier series: Setting e„ {9) = e 

/(n) := (e„ I / ) = i- T e'^V W dO, (1.1.4) 



and 



Similarly if / e (M), then 

f{t):= I e-^^'f{x)dx, (1.1.6) 



and 



Let J be an index set. We shall only need to consider the case when J 
is countable. Let {ilJa}ae 7 be a family of nonzero vectors in a Hilbert space 
7i. We say it is an orthonormal basis (ONB) if 

(V'a I V'/s) = ^a,f3 (Kronecker delta) (1.1.8) 

and if 

I]l(^a|/)l' = ll/ll' holds for all / e (1.1.9) 

aeJ 

If only (1.1.9) is assumed, but not (1.1.8), we say that {'tpa}aeJ ^ 
(normalized) tight frame. We say that it is a frame with frame constants 

0<A<B<ooii 

A Wff < 5^ |( I /)!" < S 11/11' holds for all f €H. 

Introducing the rank-one operators Qa '■= IV'a) (V'al of Dirac's terminology, 
see [BrJo02b], we see that {V'alaeJ ONB if and only if the QaS are 
projections and 

Qa = I (= the identity operator in H). (1.1.10) 

It is a (normalized) tight frame if and only if (1.1.10) holds but with no 
further restriction on the rank-one operators Qa- It is a frame with frame 
constants A and B if the operator 

S:=J2Qc, (l-l-ll) 



Februciry 1, 2008 Master Review Vol. 9in x 6in - (for Lecture Note Series, IMS, NUS) umfwaspw 



10 P.E.T. Jorgensen 

satisfies 

AI <S <BI 

in the order of hermitian operators. (We say that operators Hi = H*, 
i = l,2, satisfy ff i < ff2 if ( / | Hif ) < ( / | i?2/ ) holds for ah feH). 

Wavelets in (R) are generated by simple operations on one or more 
functions ijj in (K), the operations come in pairs, say scaling and trans- 
lation, or phase-modulation and translations. If iV <E {2, 3, . . . } we set 

V'j,fe (x) := Af^/V {N^x - k) for j, k G Z. (1.1.12) 
1.1.2. Connections to group theory 

We stress the discrete wavelet transform. But the first line in the two ta- 
bles below is the continuous one. It is the only treatment we give to the 

continuous wavelet transform, and the corresponding coherent vector de- 
compositions. But, as is stressed in [Dau92], [Kai94], and [KaLe95], the 
continuous version came first. 



Summary of and variations on the resolution of the identity operator 
1 in L"^ or in i'^, for ^ and where tpr,s {x) = r~^tp (^7^), = 

/m R (^)I' < similarly for ^ and = (a;): 



N = 2 


Overcomplete Basis 


Dual Bases 


continuous 
resolution 












discrete 
resolution 


jez feGZ 

tpj^k corresponding to 
r = 2-i,s = k2-i 


jez fcez 


N>2 


Isometrics in 


Dual Operator System in £^ 


sequence 
spaces 


N-l 
i=0 

where Sq, . . . , ^jv-i 
are adjoints to the 
quadrature mirror filter 
operators Fi, i.e.. Si = F* 


N-l 

i=0 

for a dual 
operator system 
5*0, . . . , S'at-i, 

Sq, . . . , S]SI-i 
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Consult Chapter 3 of [Kai94] for the continuous resolution, and 
Section 2.2 of [BrJo02b] for the discrete resolution. If h, k are vectors 
in a Hilbert space H, then the operator A = \h) {k\ is defined by 
the identity ( u \ Av ) = {u \ h) {k \ v) for all u,v € H. Then the 
assertions in the first table amount to: 



ff drds 

S J J ,2 

= 


\{i^r,s\f)f 

V/ e (R) 


= {f\9) 


) ( V'r,s 1 9 ) 






EE^ / 1 '^3<k) {'ii>j,k 


\9) 


jez feez 








= ll/lli^ 


V/ € (M) 


= (/lff) 


yf,geL^ (K) 


JV-l 
i=0 


\cf Vce£2 


J2{s*c\s:d) ={c 

i=0 


rf) yc,d€f 



A function ■0 satisfying the resolution identity is called a coherent vector 
in mathematical physics. The representation theory for the (ax + &)-group, 
i.e., the matrix group G = \ a £ M+, 6 G K }, serves as its underpin- 

ning. Then the tables above illustrate how the {tpj,k} wavelet system arises 
from a discretization of the following unitary representation of G: 

(f/(g.)/)(.)=«-/(^) (1.1.13) 

acting on (R). This unitary representation also explains the discretiza- 
tion step in passing from the first line to the second in the tables above. The 
functions { ipj^k \ j,k gZ} which make up a wavelet system result from the 
choice of a suitable coherent vector ij) & L? (M) , and then setting 

{x) = [u^ ,.2-. ^ {x) = 2h {2^x - k) . (1.1.14) 

Even though this representation lies at the historical origin of the subject 
of wavelets (see [DGM86]), the (aa; + 6)-group seems to be now largely 

forgotten in the next generation of the wavelet community. But Chapters 
1-3 of [Dau92] still serve as a beautiful presentation of this (now much 
ignored) side of the subject. It also serves as a link to mathematical physics 
and to classical analysis. 

Since the representation U in (1.1.13) on (M), when a unitary U is 
defined from (1.1.13) setting a = 2, 6 = 0, {Uf){x) := 2-5/ (|), leaves 
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invariant the Hardy space 

H+ = {feL^R)\ supp (/) C [0, oo) } , (1.1.15) 

formula (1.1.14) suggests that it would be simpler to look for wavelets in 
H+. After all, it is a smaller space, and it is natural to try to use the causal- 
ity features of W+ implied by the support condition in (1.1.15). Moreover, 
in the world of the Fourier transform, the two operations of the formulas 
(1.1.13) and (1.1.14) take the simpler forms 

f ^ ah^'''' f (at) and ^ ^ 2i e''^' i; {2H) . (1.1.16) 

So in the early nineties, this was an open problem in the theory, i.e., whether 
or not there are wavelets in the Hardy space; but it received a beautiful 
answer in [Aus95]. Auscher showed that there are no wavelet functions ip 
in W+ which satisfy the following mild regularity properties: 

{Ro) i) is continuous; 

{Re) for some £ e IR+, ^ {t) = O {\t\^) 

and V3 {t) =o{(l + |t|)"^"^) , t e R. 

Comparison of formulas (1.1.13) and (1.1.14) shows that The tradi- 
tional discrete wavelet transform may be viewed as the restriction to a 
subgroup H oi B. classical unitary representation of G. The unitary repre- 
sentations of G arc completely understood: the set of irreducible unitary 
representations consists of two infinite-dimensional incquivalent subrepre- 
sentations of the representation (1.1.13) on (M), together with the one- 
dimensional representations ( g J ) a**^ parameterized by fc S R. (The 
two subrepresentations of (1.1.13) are obtained by restricting to / G (K) 
with supp/ C (— c)o,0] and supp/ C [0, oo), respectively.) However, the 
subgroup H oi G has a rich variety of inequivalent infinite-dimensional 
representations that do not arise as restrictions of (1.1.13), or of any repre- 
sentation of G. The group H considered in (1.1.14) is a semidirect product 
(as is G): it is of the form 



Hn = 




a = N^, b = ^ niN\ j gZ, mG Z, 

where the summation is finite > . (1.1.17) 

(In the jargon of pure algebra, the nonabelian group H?q is the semidirect 
product of the two abelian groups Z and Z [-^j , with a naturally defined 
action of Z on Z [-^] .) 
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The papers [DaLa98], [JorOla], [BaMe99], [HLPS99], [LPTOl], and 
[BreJo91] show that it is possible to use these nonclassical representations 
of H for the construction of unexpected classes of wavelets, the wavelet sets 
being the most notable ones. Recall that a subset E cM. oi finite measure 
is a wavelet set if ip = Xe '^^ such that, for some N G Z+, N >2, the func- 
tions I N^il) {N^x — k) I j, A; e Z I form an orthonormal basis for (M). 
Until the work of Larson and others, see [DaLa98] and [HLPS99], it was 
not even clear that wavelet sets E could exist in the case N > 2. The paper 
[LPTOl] develops and extends the representation theory for the subgroups 
H]^ independently of the ambient group G and shows that each has 
continuous series of representations which account for the wavelet sets. The 
role of the representations of the groups and their generalizations for 
the study of wavelets was first stressed in [BreJo91]. 

There is a different transform which is analogous to the wavelet trans- 
form of (1.1.13)-(1.1.14), but yet different in a number of respects. It is the 
Gabor transform, and it has a history of its own. Both are special cases of 
the following construction: Let G be a nonabelian matrix group with center 
C, and let U he a unitary irreducible representation of G on the Hilbert 
space L'^ (M). When tp G L"^ (R) is given, we may define a transform 

(Tv,/) (0 := ( (0 1 / ) , for f € (R) and ^ € G/G. (LL18) 

It turns out that there are classes of matrix groups, such as the ax + 
h group, or the 3-dimcnsional group of upper triangular matrices, which 
have transforms T,p admitting effective discretizations. This means that it 
is possible to find a vector t/; e (K), and a discrete subgroup A C G/G, 
such that the restriction to A of the transform in (1.1.18) is injective 
from L"^ (R) into functions on A. 

There are many books on transform theory, and here we are only making 
the connection to wavelet theory. The book [Per86] contains much more 
detail on the group-theoretic approach to these continuous and discrete 
coherent vector transforms. 

1.1.3. Some background on matrix functions in mathematics and in 
engineering 

One of our coordinates for the landscape of multiresolution wavelets takes 
the form of a geometric index. In fact, it involves a traditional operator- 
theoretic index with values in Z. When it is identified with a winding num- 
ber or a counting of homotopy classes, it serves also as a Fredholm index of 
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an associated Toeplitz operator. An orthogonal dyadic wavelet basis has its 
wavelet function ijj satisfying the normalization ||f/'|lL2(R) = 1, i.e., '0 is a vec- 
tor of norm one in the Hilbert space (M) . In the lingo of quantum theory, 
ip is therefore a pure state, and the x-coordinate is an observable called the 
position. The integral E^p (x) = Jg_x\ip {x)f dx is the expected value of the 
position. If denotes the standard Haar function in (1.2.15), then clearly 
E^jj (x) = ^. Also note the translation formula -©^(.^fc) (x) = E^ (x) + k. 
We showed in Corollary 2.4.11 of [BrJo02b], completely generally, that the 
other orthonormal wavelets i/' have expected values in the set ^ +2. Hence, 
after V' is translated by an integer, you cannot distinguish it from the Haar 
wavelet ipH in (1.2.15) by looking only at the expected value of its position 
coordinate. The translation integer k turns out to be a winding number. 
Our result holds more generally when the definition of E^ {x) is adapted 
to a wider wavelet context, as we showed in Chapter 6 of [BrJo02b]; but in 
all cases, there is a winding number which produces the above-mentioned 
integer translate k. 

The issue of connectedness for various classes of wavelets is a general 
question which has been addressed previously in the wavelet literature; sec, 
e.g., [HLPS99], [HeWe96], [StZhOl], and [RcWe98]. Here we bring homotopy 
to bear on the question, and we identify the connected components when 
the compact support is fixed and given. We show among other things that 
for a fixed ii'i-class a homotopy may take place within a variety of wavelets 
which is specified by a slightly bigger support than the initially given one. 

An important point of our present discussion, beyond the mere fact of 
compact support, is the size of the support of the wavelets in question. 
Consider two wavelets A and B of a certain support size. Then our first 
results in this section also specify the paths C {t), if any, which connect 
A and B, and in particular the size of the support of the wavelets corre- 
sponding to C (t) . In [Br Jo02b] , we treat connectivity in the wider context 
of noncompactly supported wavelets, following at the outset [Gar99], which 
considers scale number N = 2, and wavelets V' satisfying 



Garrigos considers, for i < a < oo, the class Wa of wavelets such that 




is an orthonormal basis (ONB) for (M) . 



(1.1.19) 




(1.1.20) 
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and there is an £ = e {ip) such that 

/ ^pit)"^ (l + \tfy dt <oo, (1.1.21) 

i.e., the wavelet is supposed to have some degree of smoothness in the sense 
of Sobolev. 

We now turn to the group of functions U: T U(iV), where U (iV) 
denotes the group of all complex N-hy-N matrices. The functions will not 
be assumed continuous in general. The continuous functions will be des- 
ignated C(T,U(A^)). Each function in C(T,U(A^)) has a i^i-class, also 
called a winding number; see [BrJo02b]. The functions in C (T, U (N)) with 
finite Fourier expansion will be called Fourier polynomials, also if they are 
functions which take values in U (N). 

Proposition 1.1.3.1: Let U G C {T,l] (N)) be a Fourier polynomial, and 
assume that K\ {U) = e Z. Then U is homotopic in C (T, U (A?')) to 

V{z) = zS®{In-v) (1.1.22) 

where p is the one- dimensional projection onto the first coordinate slot in 
C^, and if U has the form 

D 

U{z)= J2 ^''^k, (1.1.23) 

k=-D 

then U may be homotopically deformed toV inC (T, U {N)) through Fourier 
polynomials of degree at most \d\ + ND. 

This proposition remains true if the word "Fourier polynomial" is re- 
placed by "polynomial" and Uk = for k = —D, —D + 1, . . . , —1. In that 
case d G Z+ and U may be homotopically deformed to V in the loop semi- 
group of polynomial unitaries in C (T, U (N)) through polynomials of degree 
at most d. 

Proof: Multiplying U by , we obtain a polynomial z^U (z) of degree 2D 
mapping T into U (N). Then Ki (z'^U) =d+ ND. By Proposition 3.3 of 
[BrJo02a], there exist d+ND onc-dimensional projections pi,p2, . . . ,Pd+ND 
in Mjv (C) and a unitary Vq € Mj^ (C) such that 

d+ND 

z'^U{z) = Vo n i^-Pi + zpi). (1.1.24) 

k=l 

(Sec § 2.2.4 for a related, but different, decomposition.) Now, deforming 
each of the pi's continuously through one-dimensional projections to the 
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projection po onto the first coordinate direction, and deforming Vq in U (iV) 
into In, we see that z^U (z) can be deformed into 

d+ND 

n {l-Po + zpo) = l-Po + z''+'"'po. (1.1.25) 

fe=i 

Thus U {z) itself is deformed into 

^-^(l-po) + ^'+('^-')^Po. (1.1.26) 

But writing (1 — po) as a sum of TV — 1 onc-dimcnsional projections 
qi, . . . , Qn-i, we have that the unitary that U {z) is deformed into is 

JV-l 

n ((1 - Qk) + z-'^q,) ■ (l + .''+(^-^)%) , (1.1.27) 

k=l 

and next deforming each of the qk in this decomposition into po, we see 
that U (z) is deformed into 

N-l 

H ((1 - po) + z-'^po) ■ (l + 0''+(^-i)^po) = (1 - Po) + zW (1.1.28) 

fc=i 

The crude estimate \d\ + ND on the degree of the Fourier polynomials 
occurring during the deformation is straightforward. 

To prove the last statement in the proposition one does not need to 
multiply U by z^ , and the proof simplifies. Note in particular that D < d 
(assuming au ^ 0). □ 

Remark 1.1.3.1: We do not know if Proposition 1.1.3.1 is true if 
C (T, U (iV)) is replaced by C (T, GL (N)). It is known from Lemma 11.2.12 
of [RLLOO] that if ^ G C(T,GL(A^)) is a polynomial of degree 1 in z, 
then A can be homotopically deformed through first-order polynomials in 
C (T, GL (N)) to a unitary of the form z — > zp + {In ^ p) for some projec- 
tion p, and hence Proposition 1.1.3.1 for C (T, GL (iV)) would follow if any 
polynomial A G C {T, GL {N)) could be factored into first-order polynomi- 
als. It is also clear, since any element A E C (T, GL (N)) can be homotopi- 
cally deformed into z'^p ® {In — p) in C (T, GL {N)), that if A is a Fourier 
polynomial, then A can be homotopically deformed into z'^p © {In — p) 
through Fourier polynomials. This follows by compactness and the Stone- 
Weierstrafi theorem (Lemma 11.2.3 of [RLLOO]). For our purposes in wavelet 
theory, though, we would need a computable upper bound for the degree 
of the Fourier polynomials. 
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For ease of reference we will now list the correspondences between the 
various objects that interest us in this case. These objects are: 

(i) matrix functions, A: T Ujv (C), satisfying the normalization 

A{1) = H, iJfc,( = -^e^2.fe;/iv, fc,/ = o,...,iV-l, (1.1.29) 

(ii) high- and low-pass wavelet filters rrii, i = 0,1, ... ,N — 1, satisfying 

mi{w)mj{w)=NSij, i,j = 0,...,N-l, (1.1.30) 

and 

mo(l) = \/iV, (1.1.31) 

(iii) scaling functions (fi together with wavelet generators ipi. 

We did not specify the continuity and regularity requirements of the func- 
tions A, rrii, (fi, tpi above. This will be done differently in different contexts 
and the classes clearly depend on these added requirements. We will now 
restrict to the case that the functions if and ijji have compact support in 
[0, oo), i.e., that A and rrii are polynomials in z. Thus z ^ A{z) is a. poly- 
nomial function with 

(Aiz))* Aiz) = l, zeT. (1.1.32) 
Scaling functions/ wavelet generators to viravelet filters {f,ip) i-^ m 

One defines a„ by 

(p{x) =VNYan>p{Nx-n), (1.1.33) 
nez 

(cf. (2.3.7)) and then mo by 

mo(^) = ^a„^", (1.1.34) 

n 

or one uses 

\^ip{Nt) =moit)ip(t) (1.1.35) 

directly. Then the high-pass filters mi,i = 1, . . . ,N—1, can be derived from 
(2.3.10) below. If we are in the generic case (2.3.6), we may also recover the 
Fourier coefficients ai*' of rrii by 

a« = (l/^/7v)(<^(. -n)\i,i{-/N)) 

= {ip{- -n)\ Uipi ) (with Vo = V>), 
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where Uipi{x) := N^^^^tpi (x/N). In particular it follows in this generic 
case that if the scaling and wavelet functions have compact support and 
the filters are Lipschitz, then the filters are Fourier polynomials. Is this true 
also in the nongeneric tight frame case? 

Now, if £) e N, define: 

• MF (D) = the set of polynomial functions in 2: G T in 

C (T, Uat (C)) of degree at most D satisfying 
(1.1.29); 

• WF (D) = the set of A^-tuples of wavelet filters 

(mo, . . . , mAT-i) such that all rrii are polynomials in 
z G T of degree at most D satisfying (1.1.30) and 
(1.1.31); 

• SF (D) = the set of A''-tuples . . . , tpN-i) of scaling 

functions/ wavelet functions with support in [0, £>]. 

The spaces MF (D), WF (D), and SF (D) may be equipped with the obvious 
topologies, coming in the first two cases from, for example, the L°°-norm 
over z, and in the last case either from the (R)-norm or, as will be more 
relevant, the tempered-distribution topology. By virtue of Proposition 3.2 
in [BrJo02a], MF (D) has the structure of a compact algebraic variety, and 
so by (2.3.4) below, WF (D) is a compact algebraic variety. It is clear from 
(2.3.4) that the map m maps MF (D) into WF {{D + 1) N - 1), and 
that m ^ A maps WF ((£) + 1) AT - 1) into MF(£)). Furthermore, it is 
clear from (1.1.33) and (2.3.10) that m {ip, ip) maps WF {{N - 1) D) into 
SF (D), and conversely {if, ^p) ^ m maps SF (D) into WF {{N - 1) D). 

Now, let a subindex denote the subsets of these various spaces such 
that the condition 

Spec(iio)nT = {1} and dim e /C^_d_j , ii (5) = = 1 (1.1.39) 

holds. It is known that the set of points such that (1.1.39) does not hold 
is a lower-dimensional subvariety of the various varieties, see Section 6 of 
[JorOlb], and hence MFo (D), WFq (D), and SFq (D) contain the generic 
points in MF (D), WF (D), and SF (£>). 

We now summarize the local connectivity results by stating the following 
theorem. The proof may be found in [BrJo02b], where this is Theorem 2.1.3. 

Theorem 1.1.3.1: Let fc G N. Equip the space SF (fcA^ + 1) of scaling 
functions/ wavelet functions with support in [0, kN + 1] with the tempered- 
distribution topology. Then SF {kN + 1) is homeomorphic to a com- 



(1.1.36) 
(1.1.37) 

(1.1.38) 
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pact algebraic variety. Furthermore, for two elements (<Po> V'o) > ('/'ii V^i) € 
SF (kN +1), the following conditions are equivalent: 

(a) The elements (i^oiV'o) and ('/'ii^i) can be connected to each other by 
a continuous path in SF {NkN + 1); 

(b) i^i((^o,V'o) = ^1(^1,^1); 

(c) The elements (tpcV'o) and {(fij-ipi) can be connected to each other by 
a continuous path in some SF (K). 

Thus, SF {kN + 1) is divided into Nk (iV — 1) + 1 components which are 
connected over SF {NkN +1). 

1.2. Motivation 

In addition to the general background material in the present section, the 
reader may find a more detailed treatment of some of the current research 
trends in wavelet analysis in the following papers: [Jor03a] (a book re- 
view), [JorOSb] (a survey), and the research papers [DuJo03], [DuJo04a], 
[DuJo04b], [DuJo04c], [Jor04a], and [Jor04b]. 

As a mathematical subject, the theory of wavelets draws on tools from 
mathematics itself, such as harmonic analysis and numerical analysis. But 
in addition there are exciting links to areas outside mathematics. The con- 
nections to electrical and computer engineering, and to image compression 
and signal processing in particular, are especially fascinating. These inter- 
connections of research disciplines may be illustrated with the two subjects 
(1) wavelets and (2) subband filtering [from signal processing]. While they 
are quite different, and have distinct and independent lives, and even have 
different aims, and different histories, they have in recent years found com- 
mon ground. It is a truly amazing success story. Advances in one area 
have helped the other: subband filters are absolutely essential in wavelet 
algorithms, and in numerical recipes used in subdivision schemes, for ex- 
ample, and especially in JPEG 2000 an important and extraordinarily 
successful image-compression code. JPEG uses nonlinear approximations 
and harmonic analysis in spaces of signals of bounded variation. Similarly, 
new wavelet approximation techniques have given rise to the kind of data- 
compression which is now used by the FBI [via a patent held by two math- 
ematicians] in digitizing fingerprints in the U.S. It is the happy marriage of 
the two disciplines, signal processing and wavelets, that enriches the union 
of the subjects, and the applications, to an extraordinary degree. While the 
Tise of high-pass and low-pass filters has a long history in signal process- 
ing, dating back more than fifty years, it is only relatively recently, say the 
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mid-1980's, that the connections to wavelets have been made. Multiresolu- 
tions from optics are the bread and butter of wavelet algorithms, and they 
in turn thrive on methods from signal processing, in the quadrature mir- 
ror filter construction, for example. The effectiveness of multircsolutions in 
data compression is related to the fact that multircsolutions are modelled 
on the familiar positional number system: the digital, or dyadic, represen- 
tation of numbers. Wavelets are created from scales of closed subspaces of 
the Hilbert space (M) with a scale of subspaces corresponding to the 
progression of bits in a number representation. While oversimplified here, 
this is the key to the use of wavelet algorithms in digital representation of 
signals and images. The digits in the classical number representation in fact 
are quite analogous to the frequency subbands that are used both in signal 
processing and in wavelets. 
The two functions 



ip{x) 



< a; < 1 
elsewhere 



and V' {x) 



1 < a; < i 

-1 i < .T < 1 

elsewhere 



(1.2.1) 



Father function Mother function 

(a) (b) 

capture in a glance the refinement identities 

ip{x) = ip{2x)+ip{2x-l) and V (a;) = V (2a;) - y (2a; - 1) . 

The two functions are clearly orthogonal in the inner product of (M), and 
the two closed subspaces Vo and Wo generated by the respective integral 
translates 



satisfy 



{ip{ - -k) ■.k&Z} and {i) {■ - k) : k ^Z} 



UVo C Vo and UWq C Vo 



(1.2.2) 
(1.2.3) 
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where U is the dyadic scahng operator Uf (x) = 2^^/^/ {x/2). The factor 
2"-^/^ is put in to make U a unitary operator in the Hilbert space L"^ (M). 
This version of Haar's system naturally invites the question of what other 
pairs of functions ip and ip with corresponding orthogonal subspaccs Vq 
and Wo there are such that the same invariance conditions (1.2.3) hold. 
The invariance conditions hold if there are coefficients Uk and bk such that 
the scaling identity 

ip{x) = J2('k'fii2x-k) (1.2.4) 
feez 

is solved by the father function, called (p, and the mother function ip is 
given by 

^{x) =^bkip{2x-k). (1.2.5) 
fcez 

A fundamental question is the converse one: Give simple conditions on two 
sequences (a^) and (bk) which guarantee the existence of (M)-solutions 
if and tp which satisfy the orthogonality relations for the translates (1.2.2). 
How do we then get an orthogonal basis from this? The identities for Haar's 
functions (p and ip of (1.2.1) (a) and (1.2.1) (b) above make it clear that the 
answer lies in a similar tiling and matching game which is implicit in the 
more general identities (1.2.4) and (1.2.5). Clearly we might ask the same 
question for other scaling numbers, for example x — > 3x or x — > 4x in place 
of X 2x. Actually a direct analogue of the visual interpretation from 
(1.2.1) makes it clear that there are no nonzero locally integrable solutions 
to the simple variants of (1.2.4), 

^{x) = ^{cp{3x)+^{3x-2)) (1.2.6) 

or 

ip{x) = 2{ip{4:x) + (p{Ax-2)). (1.2.7) 

There are nontrivial solutions to (1.2.6) and (1.2.7), to be sure, but they 
are versions of the Cantor Devil's Staircase functions, which are prototypes 

of functions which arc not locally integrable. 

Since the Haar example is based on the fitting of copies of a fixed "box" 
inside an expanded one, it would almost seem unlikely that the system 
(1.2.4)-(1.2.5) admits finite sequences (a/j) and (bk) such that the corre- 
sponding solutions If and ip arc continuous or diffcrcntiable functions of 
compact support. The discovery in the mid-1980's of compactly supported 
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and once the solutions were found, other apphcations followed at a rapid 
pace: see, for example, the ten books in Benedetto's review [BenOO]. It is the 
solution ip in (1.2.5) that the fuss is about, the mother function; the other 
one, if, the father function, is only there before the birth of the wavelet. 
The most famous of them are named after Daubechies, and look like the 
graphs in Figure 1. With the multiresolution idea, we arrive at the closed 
subspaces 

Vj := U-^Vo, j e Z, (1.2.8) 

as noted in (1.2.2)-(1.2.3), where U is some scaling operator. There are ex- 
tremely effective iterative algorithms for solving the scaling identity (1.2.4): 
see, for example. Example 2.5.3, pp. 124-125, of [BrJo02b]*, [Dau92], and 
[StNg96], and Figure 1. A key step in the algorithms involves a clever choice 
of the kind of resolution pictured in (1.2.13), but digitally encoded. The 
orthogonality relations can be encoded in the numbers (a^) and (6fe) of 
(1.2.4)-(1.2.5), and we arrive at the doubly indexed functions 

V-,- fe {x) := 2^'/2v^ {2^x -k), j,kG Z. (1.2.9) 

It is then not difficult to establish the combined orthogonality relations 

/ V'j.fe {x) i^j-,k' (x) dx = (Vj.^fc I Vj'.fe') = Sj,j'h,k' (1.2.10) 
Jm. 

plus the fact that the functions in (1.2.9) form an orthogonal basis for 
L'^ (M). This provides a painless representation of (K)-functions 

jez feez 

where the coefficients Cj^k are 

Cj,k = [ J~iMf{x) dx = <^,. I /) . (1.2.12) 

Js. 

What is more significant is that the resolution structure of closed subspaces 
of L2 (M) 

• • • C V-2 C V-i C Vo C Vi C V2 C • • • (1.2.13) 

facilitates powerful algorithms for the representation of the numbers Cj^k in 
(1.2.12). Amazingly, the two sets of numbers {ak) and (6fc) which were used 



* See an implementation of the "cascade" algorithm using Mathematica, and a "cartoon" 

of wavelets cornjjiited with it, at 

http:/ /www. math. uiowa.edu/-jorgen/wavelet_motions.pdf . 
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in (1.2.4) (1.2.5), and which produced the magic basis (1.2.9), the wavelets, 
are the same magic numbers which encode the quadrature mirror filters of 
signal processing of communications engineering. On the face of it, those 
signals from communication engineering really seem to be quite unrelated 
to the issues from wavelets — the signals are just sequences, time is discrete, 
while wavelets concern (R) and problems in mathematical analysis that 
are highly non-discrctc. Dual filters, or more generally, subband filters, were 
invented in engineering well before the wavelet craze in mathematics of 
recent decades. These dual filters in engineering have long been used in 
technology, even more generally than merely for the context of quadrature 
mirror filters (QMF's), and it turns out that other popiilar dual wavelet 
bases for (R) can be constructed from the more general filter systems; 
but the best of the wavelet bases are the ones that yield the strongest form 
of orthogonality, which is (1.2.10), and they are the ones that come from the 
QMF's. The QMF's in turn are the ones that yield perfect reconstruction of 
signals that are passed through filters of the analysis-synthesis algorithms of 
signal processing. They are also the algorithms whose iteration corresponds 
to the resolution sytems (1.2.13) from wavelet theory. 

While Fourier invented his transform for the purpose of solving the 
heat equation, i.e., the partial differential equation for heat conduction, the 
wavelet transform (1.2.11)-(1.2.12) does not diagonalize the differential op- 
erators in the same way. Its effectiveness is more at the level of computation; 
it turns integral operators into sparse matrices, i.e., matrices which have 
"many" zeros in the off-diagonal entry slots. Again, the resolution (1.2.13) 
is key to how this matrix encoding is done in practice. 



1.2.1. Some points of history 

The first wavelet was discovered by Alfred Haar long ago, but its use was 
limited since it was based on step- functions, and the step- functions jump 
from one step to the next. The implementation of Haar's wavelet in the 

approximation problem for continuous functions was therefore rather bad, 
and for differentiable functions it is atrocious, and so Haar's method was 
forgotten for many years. And yet it had in it the one idea which proved 
so powerful in the recent rebirth (since the 1980's) of wavelet analysis: the 
idea of a m,ultiresolution. You see it in its simplest form by noticing that 
a box function B of (1.2.14) may be scaled down by a half such that two 
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copies B' and B" of the smaller box then fit precisely inside B. See (1.2.14). 





B' 


B" 

















1 


2 



B X 












1 

2 


1 


X 



(1.2.15) 



This process may be continued if you scale by powers of 2 in both directions, 
i.e., by 2^ for integral fc, — oo < k < oo. So for every fc G Z, there is a finer 
resolution, and if you take an up- and a shifted mirror image down- version 
of the dyadic scaling as in (1.2.15), and allow all linear combinations, you 
will notice that arbitrary functions / on the line — oo < x < oo, with 
reasonable integrability properties, admit a representation 

/(a;) = 5^Cfc,„V(2'=a;-n), (1.2.16) 

where the summation is over all pairs of integers k,n € Z, with k rep- 
resenting scaling and n translation. The very simple idea of turning this 
construction into a multiresolution ("multi" for the variety of scales in 
(1.2.16)) leads not only to an algorithm for the analysis/synthesis problem, 

fix)^Ck,n, (1.2.17) 

in (1.2.16), but also to a construction of the single functions ip which solve 

the problem in (1.2.16), and which can be chosen differentiable, and yet 
with support contained in a fixed finite interval. These two features, the 
algorithm and the finite support (called compact support), are crucial for 
computations: Computers do algorithms, but they do not do infinite in- 
tervals well. Computers do summations and algebra well, but they do not 
do integrals and differential equations, unless the calculus problems are 
discretized and turned into algorithms. 

In the discussion to follow, the multiresolution analysis viewpoint is 
dominant, which increases the role of algorithms; for example, the so-called 
pyramid algorithm for analyzing signals, or shapes, using wavelets, is an 
outgrowth of multiresolutions. 

Returning to (1.2.14) and (1.2.15), wc sec that the scaling function cp 
itself may be expanded in the wavelet basis which is defined from tp, and 
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we arrive at the infinite series 

oo 

(p(a;) =^2-'=i/'(2-'=a;) (1.2.18) 
fe=i 

which is pointwise convergent for a; e R. (It is a special case of the expansion 
(1.2.16) when / = cp.) In view of the picture ((> ) below, (1.2.18) gives an 
alternative meaning to the traditional concept of a telescoping infinite sum. 
If, for example, < a; < 1, then the representation (1.2.18) yields (p {x) — 
1 = 5 + (i)' + ---'Whileforl<x<2,(^ (a;) =0 = -i + (i)% (!)' + •••. 
More generally, if n e N, and 2""^ < a; < 2", then 



(x) = = - 




So the function ip is itself in the space Vo C (M), and ip represents the 
initial resolution. The tail terms in (1.2.18) corresponding to 

E2-V(2-'=x)=^^(|,) (1.2.19) 

represent the coarser resolution. The finite sum 

n 
k=l 

represents the missing detail of as a "bump signal". While the sum on 
the left-hand side in (1.2.19) is infinite, i.e., the summation index k is in 
the range n < k < oo, the expression 2~"</j (2~"a;) on the right-hand side is 
merely a coarser scaled version of the original function ip from the subspace 
V C (K) which specifies the initial resolution. Infinite sums are analysis 
problems while a scale operation is a single simple algorithmic step. And 
so we have encountered a first (easy) instance of the magic of a resolution 
algorithm; i.e., an instance of a transcendental step (the analysis problem) 
which is converted into a programmable operation, here the operation of 
scaling. (Other more powerful uses of the scaling operation may be found in 
the recent book [Mey98] by Yves Meyer, especially Ch. 5, and [HwMa94].) 

The sketch below allows you to visiializc more clearly this resolution 
versus detail concept which is so central to the wavelet algorithms, also 
for general wavelets which otherwise may be computationally more difficult 
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than the Haar wavelet. 





ip{x) 








\m ^^^^^ 



The wavelet decomposition of Haar's bump function <p in (1.2.14) and (1.2.18) 
Using the sketch we see for example that the simple step function 

/ {x) = atp (x) + bip{x-l)= aX[o,i> W + &X[i,2> W (1.2.20) 




12 

has the wavelet decomposition into a sum of a coarser resolution and an 
intermediate detail as follows: 

intermediate detail coarser version 

Thus the details are measured as differences. This is a general feature that 
is valid for other functions and other wavelet resolutions. See, for instance, 
§ 2.2 below. 

1.2.2. Some early applications 

While the Haar wavelet is built from flat pieces, and the orthogonality prop- 
erties amount to a visual tiling of the graphs of the two functions (p and tp, 
this is not so for the Daubechies wavelet nor the other compactly supported 
smooth wavelets. By the Balian Low theorem [Dau92]. a time-frequency 
wavelet cannot be simultaneously localized in the two dual variables: if tp is 
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(1.2.23) 



a time-frequency Gabor wavelet, then the two quantities \xip {x)f dx and 

ttj' (t) (it cannot both be finite. Since (^jj^ (t) = iti' (t), this amounts 
to poor differentiabihty properties of well-localized Gabor wavelets, i.e., 
wavelets built using the two operations translation and frequency modula- 
tion over a lattice. 

Tint with the multircsohition viewpoint, wc can understand the first of 
Daubechies's scaling functions as a one-sided diffcrcntiablc solution tp to 

(fi (x) = ho(p {2x) + hup {2x + h2(p {2x -2) + h^ip {2x - 3) , (1.2.22) 

where the four real coefficients satisfy 

ho + hi+h2 + h3 = 2,' 

hs- h2 + hi - ho = 0, 

hs - 2/i2 + Shi - 4/io = 0, 

hihs + hoh2 = 0. 

The system (1.2.23) is easily solved: 

4/io = 1 + VS, 4/i2 = 3 - 

4/ii = 3 + VS, 4/i3 = 1 - 

and Daubcchies showed that (1.2.22) has a solution ip which is supported 
in the interval [0, 3] , is one-sided diffcrcntiablc, and satisfies the conditions 

ip{x) dx = 1, / V {x) dx = 0, and / xip (x) dx = 0. (1.2.25) 

Js. Jr Jr 

The first applications served as motivating ideas as well: optics, seis- 
mic measurements, dynamics, turbulence, data compression; see the book 
[KaLe95] Actually, it is two books: the first one (primarily by Kahane) is 
classical Fourier analysis, and the second one (primarily by P.-G. Lemarie- 
Rieusset) is the wavelet book. It will help you, among other things, to get 
a better feel for the French connection, the Belgian connection, and the 
diverse and early impulses from applications in the subject. Enjoy! 

For a list of more recent applications we recommend [MeyOO]. 



(1.2.24) 



2. Signal processing 

If we idealize and view time as discrete, a copy of Z, then a signal is a se- 
quence (^n)„£z of numbers. A filter is an operator which calculates weighted 
averages 

(U^l^afcCn-fc- (2-1) 

feez 
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Fig. 2. Perfect reconstruction of signals 



But working instead with functions of z G T, this is muhiphcation, / (z) i— > 
m (z) / (z), where m (z) = J2kez ^k^'^ ^^'^ / (^) ~ Sfcez ^kz'' are the usual 
Fourier representation of the corresponding generating functions. Similarly, 
down-sampling @ and up-sampling (m) as operators on sequences take the 
form 



1 

N 



E 



and 



(2.2) 



(2.3) 



Since the operators (m) and (m) arc clearly dual to one another on the 
Hilbert space £^ (Z) of sequences (i.e., time-signals), we get the correspond- 
ing duality for (T), i.e., 



/ / (z^) g (z) d„ {z) = j f{z)^Y. 9 i'^) dl^ {z) , (2. 



4) 



where /U denotes the normalized Haar measure on T, or equivalently the 
following identity for 27r-periodic functions: 



'9 + k-2TT 



de. 



(2.5) 



Quadrature mirror filters with N frequency subbands mo, mi, . . . ,mN-i 
give perfect reconstruction when signals are analyzed into subbands and 
then reconstructed via the Tip-sampling and corresponding dual filters. In 
engineering formalism this is expressed in the diagram in Fig. 2, for N = 2, 
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and mo, resp. mi, are called low-pass, resp. high-pass, filters. In operator 
language, this takes the form 

F*Fo + F*F, = /, 

where Fq and Fi arc the operators in Fig. 2, with dual operators Fq and 
F^ . The quadrature conditions may be expressed as 

FoFo* = F^F* = I (2.6) 

and 

FoF^ = FiF^ = 0. (2.7) 

In operator theory there is tradition for working instead with the operators 

Sj :~ F* . When viewed as operators on (T) they are therefore isometries 
with orthogonal ranges, and they satisfy 

Af-l 

E ^^^1 = ^ (2-8) 

j=0 

with / now representing the identity operator acting on (T). The rela- 
tions on the S'j-operators are known as the Cuntz relations because of their 
use in C* -algebra theory; see [Cun77]. In the present application they take 
the form 

iS,f) (z) = m, {z) f (z") , / e (T) , (2.9) 

and 

= ^ E ^M^/H' (2-10) 

—z 

and the Cuntz relations are equivalent to the conditions 

E \m,{w)f = N (2.11) 

and 

E mj (w) mk (w) = for all z e T and j 7^ k. (2.12) 

The last conditions are known in engineering as the quadrature conditions 
for the subband filters m,Q, mi. .... m,iv_i, with mo denoting the low-pass fil- 
ter. The low-pass and band-pass conditions on the functions rrij are perhaps 
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more familiar in the additive notation given by the substitution z := e~'^. 
Then the functions mj are viewed as 27r-periodic, and 



while 




with both of the indices j, k ranging over 0,1,..., A'' — 1. 
2.1. Filters in communications engineering 

The coefficients of the functions rrij {■) are called impulse response coef- 
ficients in communications engineering, and when used in wavelets and 
in subdivision algorithms, they are called masking coefficients. In the finite 
case, the ruj {■) 's are also called FIR for finite impulse response. The model 
illustrated in Fig. 2 is used in filter design in either hardware or software: 

|1] Try filters mo, mi in Fig. 2, and approximate the output to the input; 
|2| Choose a specific structure in which the filter will be realized and then 

quantize the coefficients, length and numerical values; 
|3] Verify by simulation that the resulting design meets given performance 

specifications. 

Once filters are constructed, we saw that they are also providing us 
with wavelet algorithms. When the steps of Fig. 2 are iterated, we arrive at 
wavelet subdivision algorithms. Relative to a given resolution (pictured as a 
closed subspace Vi, say, in (K)), signals, i.e., functions in (M), decom- 
pose into coarser ones and intermediate details. Relative to the subspaces 
Wo and Vi , this amounts to 

Vi = Vo + Wo. (2.1.1) 

T T T 

given coarser intermediate 

resolution resolution detail 

Ideally, we wish the decomposition in (2.1.1) to be orthogonal in the sense 
that 



( / I c, ) = for aU / e Vo and all g S Wo. (2.1.2) 
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Since the subdivisions involve translations by discrete steps, we specialize 
the resolution such that both of the spaces Vo and Wo are invariant under 
translations by points in Z, i.e., such that 

T:/^/(.-l) (2.1.3) 

leaves both of the subspaces Vo and Wo invariant. The multiresolution anal- 
ysis case (MRA) corresponds to the setup whcni Vo is singly generated, i.e., 
there is a function (/? G Vo such that the closed linear span of 

Tn^{-)^'f{--n), neZ, (2.1.4) 

is all of Vo- If = 2, then there is then also stt/j E Wq such that the closed 
linear span of { V' ( • — n) : n G Z } is all of Wo- If iV > 2, we may need func- 
tions ipi,. . . , tl>N-i in Wo such that {tpi{- — n) : i = 1, . . . ,N — 1, neZ} 
has a closed span equal to Wq. 



2.2. Algorithms for signals and for wavelets 

The pyramid algorithm and the Cuntz relations. Since the two 
Hilbert spaces (T) and i"^ (Z) are isomorphic via the Fourier series 

representation, it follows that the system {Sij^^Q is equivalent to a sys- 
tem {Si}^^Q acting on €^(Z). Specifically, {Sif) — Sif, « = 0, 1, where 
/(n) := /,j.2~"/(2;) dii{z). For c := {cn)nei ^ (^)' ^"^^ functions / on 
M, set 

/_i(x) := ([//) (x)=2-5/(|),and 
(c*/) (x) •.= ^Cnf {x - n) . 

n^TL 

For the present, let {mj}|^Q be the low-pass and high-pass wavelet filters, 
and let V' be the corresponding scaling function, resp., wavelet function, 
also called father function, resp., mother function. Now introduce the corre- 
sponding operators Si and their cousins 5,. The adjoints 5* are also called 
fxlters. 
Then 

c * = ((5*0*0) * ^) + (('^i'') * ^^^^ • 

V ' V ' 

coarser resolution detail 

Define W-.l'^^i^hy 

W (c) (a;) = {c*<fi) (x) = ^ c„(p {x - n) . (2.2.2) 

ra6Z 
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Then W maps i'^ isometrically onto Vo in the orthogonal case and 

WSo = UW. 

Further 

WSoSqC = [SqC * v?) ^ . 
Embedding into P ® as f' 0, extend to ^ f. putting 

W {c® d) = c* ip + d*ij:. 
Then the extended W maps £^ ® £^ isometrically onto U~^Vo and 
W {Soc + Sidj = UW (c ® d) 

for all c.d E P , where the left W is the one from (2.2.2) and the right is 
the extension of W to £'^®i'^. 

At this point you can use 1^2 = SqSq + SiS^ to show (2.2.1). Note that 
if Co = a and ci = b and Cj = for other i, the formula (2.2.1) reduces to 
(1.2.21). 

The subdivision relations (2.2.1) are equivalent to the system 

V2ip {2x) = a2kV {x + k) + Y hki^ [x + k), (2.2.3) 

fcez feez 

V2ip {2x - 1) = X! a2fe+i¥' {x + k) + Y^ WiV- (a; + fc) , (2.2.4) 
feez feez 

where the coefficients a„, hn are those of the quantum wavelet algorithm, 
i.e., the coefficients in the "large" unitary matrix (2.2.5). Thus the quan- 
tum algorithm does the wavelet decomposition within a fixed resolution 
subspace. 

The scaling function Lp defines a resolution subspace Vo C (M). Then 
(2.2.1), or equivalently (2.2.3)-(2.2.4), represents the orthogonal decompo- 
sition of functions in Vq into an orthogonal sum of a function with coarser 
resolution and a function in the intermediate detail subspace. 

Let 777,0. iTii 1^(- a dyadic wavelet filter, and let T 9 z i-^ A (z) S U2 (C) 
be the corresponding matrix function, Ajj (z) = 5 X)^2__j, u'~''mj (w). If 
the low-pass filter mo {z) = ao + a\z -|- • • • -|- a2n+iz'^"'^^ 1 then a choice 
for mi{z) = Yl^j^^bkZ^ is hk = {—1)'^ d2n+i-k- We then have A{z) = 

Efr-o ^kz'' where Ak = f ""^ ^'''+'1 , and the following 2"+^ x 2"+^ scalar 

\b2k 02k+lJ 
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matrix can be checked to be unitary: 



ai 
h 


Ai 


A2 




An^l 
















ao 
bo 






A„ 


Ai 




An-2 


An-1 


An 























Aa 




An-3 


An-2 


An-1 


An 






























































































0'2n+l 
&2n+l 














Ao 


Ai 


A2 




An-1 


l2n 
b2n 


<l2n-l 
b2n-l 


An 














Ao 


Ai 




An-2 


02n-2 
b2n-2 


0-271-3 
b2n-3 


An^l 


An 














Ao 




An-3 


<l2n-4 
&2n-4 


























as 
b3 


A2 


A3 




An 














Ao 


ffl2 
62 



(2.2.5) 



Except for the scalar entries in the two extreme left and right columns, 
all the other entries of the big combined matrix Ua are taken from the 
cyclic arrangements of the 2x2 matrices of coefficients Ao,Ai,...,An in 
the expansion of A{z). For the case of n = 1 this amounts to the simple 
8x8 wavelet matrix 



^0 



h 


_0_ 







0-3 

\b3 



Ai 



Ac 



Ai 



Ao 



Ai 



Ao 



ao\ 
bo 



a2 
62/ 



(2.2.6) 



Ai, 
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which is the one tliat produces the sequence of quantum gates. The quantum 
algorithm of a wavelet filter is thus represented by a 2"+^ x 2"+-^ unitary 
matrix Ua acting on the quantum qubit register C ® ■ y ® C = C^("+^', 

n+2 times 

i.e., it acts on a configuration of n + 2 qubits. The realization of a wavelet 
algorithm in the quantum realm thus amounts to spelling out the steps in 
factoring Ua into a product of qubit gates. By Shor's theorem, we know 
that this can be done, and Ua may be built out of one-qubit gates and 
CNOT gates following the ideas sketched above. The reader may find more 
discussion of the matrix Ua in Section 3 of [Pre02]. 

The generalization of classical and quantum wavelet resolution algo- 
rithms from iV = 2 to > 2 is immediate: Then m.i (z) ~ X^feez '^fe'*-^'^' 

{S,f){z)=mi{z)f{z^), i = 0,...,N-l, (2.2.7) 

and the transformation rules 

^Nk+i = i = 0,l,...,N-l, (2.2.8) 

permute the set of ONB's in (Z) and define a unitary commuting with the 
A/'-shift. Hence, the standard formulas from [Wic93], [Kla99], and [FiWi99] 
for the quantum computing algorithm naturally generalize to the case N > 
2 via (2.2.8). Instead of /e-registers (g) ■ ■ ■ (g) C^^ = C'^" over C^ we will 

k times 

now have to work rather with ® • • • = . 

k times 

The use of the algorithmic relations in engineering and operator algebra 
theory predates their more recent use in wavelet theory and wavepacket 
analysis. 

2.2.1. Pyramid algorithms 

For N > 2, the algorithm of the previous section takes the following form. 

The pyramid algorithm and the Cuntz relations revisited. By 

Fourier equivalence of (T) and (Z) via the Fourier series, it fol- 
lows that the system {<S'i}^o^ is equivalent to a system {Si}^^^ acting 
on ^2(Z). Specifically, (sj)" = Sj, i = 0,...,N -1, where /(n) := 
Jj z~"f (z) dfjL {z). For c := (c„)„g2 ^ (^)' ^"^^ functions / on M, set 
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and 



(c*/) (x) :=^c„f{x-n). 



Let {wij^o^ low-pass and high-pass wavelet filters, and let (f, 
ipi, . . . ,ippf-i be the corresponding scaling function, resp., wavelet func- 
tions. Now introduce the corresponding operators Si, and their cousins S^. 
The adjoints are also called filters. 
Then 

c*ip = ((s-t^c) *ip^ {{^i^) * ^i) aU cef (Z) . (2.2.9) 

^ .. ^ i=i 



coarser resolution 



detail 

The scaling function (p defines a resolution subspace Vo C (R). For 
the case N > 2: 

Discrete vs. continuous wavelets, i.e., £^ vs. (E): 



{0} 



V2 



Vi 



Vo 



finer scales 



W3 



W2 



Wi 



• • • rest of i 



W 



{0} 



u 


u 


U 


I 

So 


1 

So 


So 



Sic 



SoC 



W 



1=1 



sie 



More refined pyramid algorithms yield wavelet packets as follows. 

The Haar wavelet is supported in [0, 1], and if j G Z+ and k G Z, then 
the modified function x ijj (2^x — k) is supported in the smaller interval 
^ < X < ^ti. When j is fixed, these intervals are contained in [0, 1] for 
fee {0,1, ...,2-' — 1}. This is not the case for the other wavelet functions. 
For one thing, the non-Haar wavelets if) have support intervals of length 
more than one, and this forces periodicity considerations; see [CDV93]. For 
this reason, Coifman and Wickerhauser [CoWi93] invented the concept of 
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wavelet packets. They are built from functions with prescribed smoothness, 
and yet they have localization properties that rival those of the (discontin- 
uous) Haar wavelet. 

There are powerful but nontrivial theorems on restriction algorithms 
for wavelets tpj^k (x) = (2-'x - k) from (r) to (0, 1). We refer the 
reader to [CDV93] and [MiXu94] for the details of this construction. The 
underlying idea of Alfred Haar has found a recent renaissance in the work 
of Wickerhauscr [Wic93] on wavelet packets. The idea there, which is also 
motivated by the Walsh function algorithm, is to replace the refinement 
equation (1.1.33) by a related recursive system as follows: Let mo (z) = 
Y.k^kz'', mi{z) = J2k^kz'', for example bk = {-!)'' ai^k, fc € Z, be 
a given low-pass/high-pass system, N = 2. Then consider the following 
refinement system on M: 

W2n{x) = V2Y,akWn{2x-k), (2.2.10) 

W2n+i{x) = V2Y,bkWn{1x-k). (2.2.11) 

feez 

Clearly the function Wq can be identified with the traditional scaling func- 
tion Lp of (2.3.7). A theorem of Coifman and Wickerhauscr (Theorem 8.1, 
[CoWi93]) states that if P is a partition of {0,1,2,...} into subsets of the 
form 

h,n = {2^=71, 2'=n + l,...,2'= (n+l)-l}, 

then the function system 

{2*W„(2'=a;-;) | /fc,„ e « € z} 

is an orthonormal basis for (M) . Although it is not spelled out in 
[CoWi93], this construction of bases in L"^ (R) divides itself into the two 

cases, the true orthonormal basis (ONB), and the weaker property of form- 
ing a function system which is only a tight frame. As in the wavelet case, 
to get the ■P-system to really be an ONB for LP' (R), we must assume the 
transfer operator J^|„iq|2 to have Perron-Frobenius spectrum on C (T). This 
means that the intersection of the point spectrum of R\mo\'^ with T is the 
singleton A = 1, and that dimker((l — R\mo\'^)\c (j)) = 1- 

2.2.2. Subdivision algorithms 

The algorithms for wavelets and wavelet packets involve the pyramid idea as 
well as subdivision. Each subdivision produces a multiplication of subdivi- 
sion points. If the scaling is by N , then j subdivisions multiply the number 
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of subdivision points by . If the scaling is by a x d integral matrix 
N, then the multiplicative factor is IdetNl"* in the number of subdivision 
points placed in R'^. 

In the discussion below, we restrict attention to = 1, but the conclu- 
sions hold with only minor modification in the general case of d > 1 and 
matrix scaling. 

If is a continuous function on T, the transfer operator or kneading 
operator Rw 

Rw^ {z) = ^ J2W (w) ^ {w) = S*oW^ (z) , (2.2.12) 
with the alias 

(i?M./)„ = 5^CjVn-fe/fe (2.2.13) 

k 

in the Fourier transformed space, has an adjoint which is the subdivision 
operator or chopping operator 

(i?U)W=W>)^(z~) (2.2.14) 
on functions ^ on T, with the alias 

{Rwf)n = J2^^^f'^ (2.2.15) 

k 

on sequences. 

We will analyze the duality between Rw and R^ and their spectra. 
Specializing to W — |?Tio|^, we note that R\y is then the transfer operator 
of orthogonal type wavelets. In the following, W is assumed only to satisfy 
W e Lip 1 (T) and W >0. Other conditions are discussed in [BrJo02b]. 

In the engineering terminology of § 2.2, the operation (2.2.13) is com- 
posed of a local filter with the numbers Ck as coefficients, followed by the 
down-sampling (m), while (2.2.15) is composed of up-sampling (m), followed 
by an application of a dual filter. In signal processing, (m) is referred to as 
"decimation" even if N is not 10. 

The operator S (= R^) is called the subdivision operator, or the wood- 
cutter operator, because of its use in computer graphics. Iterations of S will 
generate a shape which (in the case of one real dimension) takes the form of 
the graph of a function / on R. If ^ G (Z) is given, and if the differences 



Dn{i) = f(^)-{S-0(i), ieZ, (2.2.16) 
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are small, for example if 

lim sup|£)„(i)| =0, (2.2.17) 

then we say that ^ represents control points, or a control polygon, and the 
function / is the limit of the subdivision scheme. 

It follows that the subdivision operator S on the sequence spaces, espe- 
cially on £°° (Z), governs pointwise approximation to refinable limit func- 
tions. The dual version of S, i.e., R = S* (= the transfer operator) governs 
the corresponding mean approximation problem, i.e., approximation rela- 
tive to the (IR)-norm. 

In Scholium 4.1.2 of [BrJo02b], we consider the eigenvalue problem 

= A^, A e C, (2.2.18) 

and ^ 7^ in some suitably defined space of sequences. The formula (2.2.16) 
for the limit of a given subdivision scheme S makes it clear that the case 
(2.2.18) must be excluded. For if (2.2.18) holds, for some A G C, and some 
sequence ^ of control points, then there is not a corresponding regular 
function / on M with its values given on the finer grids 2~"Z, n = 1,2,..., 

by 

(i2-") « (5"0 (i) = (i), i e Z. (2.2.19) 

We show in Example 4.1.3 of [BrJo02b] that there are no such control points 
C in ^2 (Z) \ {0}. Hence the stability of the algorithm! 

2.2.3. Wavelet packet algorithms 

The main difference between the algorithms of wavelets and those of wavelet 
packets is that for the wavelets the path in the pyramid is to one side only: 
a given resolution is split into a coarser one and the intermediate detail. 

The intermediate detail may further be broken down into frequency bands. 
With the operators Sjf {z) = m.j (z) f {z'^) acting on (T), the coarser 
subspace after j steps is modelled on S^L"^ (T), and the projection onto 
this subspace is SqSq'' where So is the isometry of (T) = Vo defined by 
the low-pass filter toq. But in the construction of the wavelet packet, the 
subspace resulting by running the algorithm j times is Si^ Si^ ■ ■ ■ Si- (T) , 
and the projection onto this subspace is 
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lin gZ.^., the wavelet function Wn is computed from the iteration ii,. . . ,ij 
corresponding to the representation 



n 



ii + i2N + izN"^ + ■■■ + ijN^^\ 



where ii, . . . ,ij € {0,1,...,7V — 1} are unique from the Euchdean algo- 
rithm. 

2.2.4. Lifting algorithms: Sweldens and more 

The discussion centers around the matrix functions A : T — !■ GL2 (C) . 

The case detA = 1. Recall that we call a finite sum J2k=-no ■^kz'', 
no,n\ > 0, a Fourier polynomial both if the coefficients Ak are numbers, 
and if they are matrices. The matrix- valued Fourier polynomials T 9 z i-^ 
A {z) € Ma (C) such that det A{z) = l form a subgroup of C (T, GL2 (C)) 
which we denote SC^- 

For every A {z) in SC2 there are m G Z+, K G C\{0}, and scalar- valued 
Fourier polynomials ui{z),. . . , Um{z), h{z),. . . , lm{z) such that 



A{z) 



(2.2.20) 



"if \ / 1 0\ /I ui (z)^ 
K-')\h{z)l)\0 1 

1 0\ flu2iz)\ f 1 0\ flumiz)^ 

See [DaSw98]. This is the first step in the Daubechies-Sweldens lifting algo- 
rithm for the; discrete wavelet transform. Thus the case det {A {z)) = 1 gives 
a constriietive lifting algorithm for wavelets, and such an algorithm has not 
been established in the C (T, GL2 (C)) case. The decomposition could also 
be compared with Proposition 3.3 of [BrJo02a], which was mentioned in 
connection with the proof of (1.1.24). 

Recall the correspondence between matrix functions and wavelet filters: 
If A : T — > GL2 (C) is a matrix function, then the corresponding dyadic 
wavelet filters are 

1 

It follows that the two matrix functions A and B satisfy 

-GO- 

for some I in the ring T of Fourier polynomials if and only if niQ^^ = m,(f^ 
and m^^^ (z) = m^' (z) + I (z^) TOq^^ (z). 
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Similarly note that the two matrix functions A and B satisiy 




for some u G if and only if m{ = m\ and mg (z) = TOq (z) + 



Remark. The conclusion is that the wavelet algorithm for a general wavelet 
filter corresponding to a matrix function, say A, may be broken down in a 
sequence of zig-zag steps acting alternately on the high-pass and the low- 
pass signal components. 

2.3. Factorization theorems for matrix functions 

We mentioned that for matrix functions corresponding to finite impulse 
response (FIR) filters which are unitary, we need only the constant matrix 
(which is choscin such as to achieve the high-pass and low-pass conditions) 
and factors of the form 



where P is a rank-one projection in C and N is the scaling number of the 
subdivision. 

Unfortunately, no such factorization theorem is available for the non- 
unitary FIR filters. But the matrix functions take values in the non-singular 
complex N X N matrices. The Sweldens-Daubechies factorization and the 
lifting algorithm serve as a substitute. There are still the general non- 
unimodular FIR-matrix functions where factorizations are so far a bit of 
a mystery. The matrix functions are called polyphase matrices in the engi- 
neering literature. The following summary serves as a classification theorem 
for the orthogonal wavelets of compact support: the wavelets correspond to 
FIR polyphase matrices which are unitary. 

In summary, an algorithm to construct all the wavelet functions tp of 
scale 2 with support in [0, 2k + 1] can be established as follows: 

|1] Pick k one-dimensional orthogonal projections Qi,...,Qk in M2 (C) 
and define the unitary-valued matrix function A (z) on T by 



A{z) = V{l-Qi + zQi) {I-Q2 + ZQ2) ■■■{1-Qk + zQk) , (2.3.1) 



u 




where 




(2.3.2) 
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Then each Qj has the form 



where € [0, 1] and Oj € [0,27r). (See Proposition 3.3 of [BrJo02a].) 
|2] Define the filters mo {z) and mi{z) by 

N-l 

mi{z)=Y,^'^iji^'')^ i,i = 0,...,7V-l, (2.3.4) 

with A'' = 2. 
|3] Define by 

If the condition 

PER (|(^|^) (i) := ^ |(^ (i + 27rn)|^ = 1 (2.3.6) 
nez 

fails, then the algorithm stops. 
|4] If the condition (2.3.6) holds, one may alternatively define (p by the 
cascade algorithm 

ip{x) = y/N'^a„(fi{Nx -n) , (2.3.7) 

, , / 1, < a; < 1, „ 

^^") = lo,.eiR\[o,i), 

Ma:tp< — > ViV ^ a„ V (A^a; - n) . (2.3.9) 

n 

[5] The wavelet function ip is then defined by 

il)i{x) = y/N^a';i'>ip{Nx-n), (2.3.10) 

where are the Fourier coefficients of m^, 

mi(^) = ^a«2", (2.3.11) 

and a; = e~'*; this is the most general wavelet function with support in 

[0,2fc+ 1]. 

|6] All other wavelet functions with compact support can be obtained from 
the ones in |5] by integer translation. 



Februciry 1, 2008 Master Review Vol. 9in x 6in - (for Lecture Note Series, IMS, NUS) umfwaspw 



Unitary Matrix Functions, Algorithms, Wavelets 43 

2.3.1. The case of polynomial functions [the polyphase matrix, joint 

work with Ola Braiteli] 

One problem occurring in the biorthogonal context which does not have 
an analogue in the orthogonal setting stems from the fact that the duality 
relations 

^ rui {w) rhj (w) = N5ij hv i,j = 0, . . . ,N - 1 (2.3.12) 

do not give any absolute restrictions on the size of rUi and rhj, e.g., a bound 
on the inner product of two vectors in does not give a bound on the 
size of the vectors if they arc not equal. This is reflected in the bi-Cuntz 
relations defined by m^, rhi. Let us now define 

(Sif) {z) = rrii {z) f (z^) , {SJ) (z) = rfn {z) f (^^) (2.3.13) 

for z € T, / e Z/^ (T). Instead of the usual Cuntz relations, the Si, Si now 
satisfy 

S:Sj=dijl, (2.3.14) 
Y^SiSt = l. (2.3.15) 

i 

If A, A G C (T, GLff (C)) are the matrix- valued functions associated to m, 
and rhi by 

m{z) = A{z^)v{z), m{z) = A{z^)v{z), (2.3.16) 
we compute 

S:S, = (AA*)^^, (2.3.17) 

in the sense that S*Sj is contained in the commutative algebra of mul- 
tiplication operators on L'^ {T) defined by C(T), and {AA*)j ^ e C(T). 
Correspondingly, 

StSj = {AA*).. (2.3.18) 

so all the operators S*Sj, S*Sj are contained in the abelian algebra C (T). 
We may introduce operators S, S from 

(T)^ = (T) . . . ^2 ^T) (2.3.19) 

N-l 

into L2 (T) by 

S = {So, Si, . . . , Sn-i) , S = {So, . . . , Sn-i) (2.3.20) 
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and then S* maps (T) into (2.3.19), etc., and the relations 
(2.3.14)-(2.3.18) take the form 

S*S =1, where 1 is the identity in Mn (C) (g) C (T) , /o 3 21^ 
SS* = 1, where 1 is the identity in C (T) , ( • • ) 

C q*q _ A A* 

These relations say that all combinations of products of S and S* with S 
and S* lie in the algebra Mjv (C) (g) C (T). But in addition A and A are 
matrix- valued functions on T, so 

AA*AA* = AA* = 1 = AA*AA* (2.3.23) 

and hence 

S*S = (5*5') ^ (2.3.24) 

and all the matrix-valued functions commute. 

This discussion can be summarized by saying that the bi-Cuntz relations 
are much less rigid than the original Cuntz relations, i.e.: 

Scholium 2.3.1.1: Given any bijcctivc operator 5* from L"^ (T)^ into 
(T) one may define S = (S*) and the bi-Cuntz relations (2.3.21) are 
satisfied. If, more specifically, S is given by (2.3.20) and (2.3.13), then oper- 
ators 5^0, . . . , Sn-1 exist such that the bi-Cuntz relations (2.3.14)-(2.3.15) 
are satisfied if and only if the operator A E Mm (C) ® C (T) defined by 
(2.3.16) is invertible, in which case one must use A = {A*)~ , (2.3.16), and 
(2.3.13) to define So,..., Sn-i- 

Let us now connect the filters to the wavelets. We have already defined 
the scaling functions if, (p and wavelet functions V',;- / = 1, . . . , A''. The 
expansions for 93 and (p converge uniformly on compacts, thus (p and Cp are 
continuous functions on M. To decide that these functions are in L? (M) one 
again forms 

{t) = PER (jpf) (t) = Y.\p{t- 2nn)f (2.3.25) 

nez 

and similarly, and one deduces again from the nonlinear intertwining 
relation 

J?'(p(V'i,V'2))=p(M^„V'i,^CV'2), ken (2.3.26) 

that 

Rmo iff) = fv ^rho if-fi) = fv- (2.3.27) 
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2.3.2. General results in mathematics on matrix functions 

In the standard case of the good old orthogonal wavelets in (R) of 
subbands, we will look for functions ipi, . . . , tpN-i in L'^ (R) such that, if k 
and n run independently over all the integers Z, i.e., — oo < fc, n < oo, then 
the countably infinite system of functions 

{ AT'^/Vi {N'^x -n)\i = l,...,N-l, k,n€Z} (2.3.28) 

is an orthonormal basis in the Hilbert space (K). The second half of 
the word "orthonormal" refers to the restricting requirement that all the 
functions ijji, . . . , ipN-i satisfy 

/ IVi {x)f dx = 1, (2.3.29) 
Jr 

or stated more briefly, 

ll^illL^(R) = 1; (2-3.30) 

or yet more briefly, 

= 1. (2.3.31) 

From familiar properties of the Lebesgue measure on R, it then follows that 
all the functions 

1p^^k,n{x) := N''/'^iP,{N''x-n), l<i<N,k,nGZ, (2.3.32) 
satisfy the normalization, i.e., that 

||V'i,fc,n|| = 1 foralli,A;,n 
The functions (2.3.32) are said to be orthogonal if 

/ i^i,k,n (x) i)i',k',n' {x) dx = Q (2.3.34) 

jR 

whenever {i,k,n) ^ {i',k',n'). We say that the two triple indices are dif- 
ferent if i ^ i' OT k ^ k' or n ^ n' . If, for example, i = i' and k = k', 
then when the same function is translated by different amounts n and n', 
the two resulting functions are required to be orthogonal. It is an elemen- 
tary geometric fact from the theory of Hilbert space that if the functions 
in (2.3.32) form an orthonormal basis, then for every function f & L? (R), 
i.e., every measurable function / on M such that 

11/11' = / |/(a;)|'da;<oo, (2.3.35) 



(2.3.33) 
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we have the identity 

r 2 

\\ff=J2 i'i,k,n{x)f{x)dx , (2.3.36) 

where the triple summation in (2.3.36) is over all configurations 1 < « < iV, 
k,n gZ. It is convenient to rewrite (2.3.36) in the following more compact 
form: 

ll/f =EKV'a,n|/)|'. (2.3.37) 

Surprisingly, it turns out that (2.3.37) may hold even if the functions ^i^k,n 
of (2.3.32) do not form an orthonormal basis. It may happen that one 
of the initial fimctions tpi, . . . , or tpN-i satisfies Wi'iW < 1, and yet that 
(2.3.37) holds for all / e (M). These more general systems are still called 
wavelets, but since they are special, they are referred to as tight frames, as 
opposed to orthonormal bases. In either case, we will talk about a wavelet 
expansion of the form 

f{x) = Y,{ ^iAn I / ) i^i,k,n {x) . (2.3.38) 

i,k,n 

It follows that the sum on the right-hand side in (2.3.38) converges in the 
norm of (M) for all functions / in (M) if (2.3.37) holds. 

But there is a yet more general form of wavelets, called biorthogonal. The 
conditions on the functions ipi, . . . , ipN-i are then much less restrictive than 
the orthogonality axioms. Hence these wavelets are more flexible and adapt 
better to a variety of applications, for example, to data compression, or to 
computer graphics. But the biorthogonality conditions are also a little more 
technical to state. We say that some given functions tpi, i = 1, . . . , A^ — 1, in 

(M) are part of a biorthogonal wavelet system if there is a second system 
of fimctions ijj,, i = 1, . . . , N - 1, in L'^ (R), such that every f £ L'^ (K) 
admits a representation 

fi^) = ^{ '^hk,n I / ) '^i,k,n {x)=^{ ■4'i,k,n \ f ) V'i,fe,n {x) , (2.3.39) 
z,fc,n i,k,n 

and 

^i,k,n {x) = N'^'^i {N'^x - n) . (2.3.40) 

In the standard normalized case where ('^j | V'j ) = 1, then you will notice 
that condition (2.3.37) turns into 

11/11' = i^i^nlf) { i>i,k,n I / ) (2.3.41) 
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for all f eL^ (M). 

The orthogonal wavelets correspond to matrix functions T ^ Uat (C), 
while the wider class of biorthogonal wavelets corresponds to the much big- 
ger group of matrix functions T GLj\t (C), via the associated wavelet 
filters. You may ask, why bother with the more technical-looking biorthog- 
onal systems? It turns out that they are forced on us by the engineers. They 
tell us that the real world is not nearly as orthogonal as the mathematicians 
would like to make it out to be. There is a paucity of symmetric orthogonal 
wavelets, and symmetry ("linear phase") is prized by engineers and work- 
ers in image processing, where the more general wavelet families and their 
duality play a crucial role. Now what if we could change the biorthogonal 
wavelets into the orthogonal ones, and still keep the essential spectral prop- 
erties intact? Then everyone will be happy. This last chapter shows that it 
is possible, and even in a fairly algorithmic fashion, one that is amenable 
to computations. 

Wavelet filters may be understood as matrix functions, i.e., functions 
from the one-torus T c C into some group of invertible matrices. If the scale 
number is N, then there are three such matrix groups which are especially 
relevant for wavelet analysis: 



Uat (C): all 
unitary N x N 
complex 
matrices 



C 



GLat (C): all 

invertible 

N X N complex 

matrices 



D 



SLn (C): all 
N X N complex 
matrices A with 
detA= 1. 



It is possible to reduce some questions in the GLjv case to better understood 
results for XJn (C); see Chapter 6 of [BrJo02b]. The SL2 case is especially 
interesting in view of Daubechies-Sweldens lifting for dyadic wavelets; see 
S 2.2.4. 



2.3.3. Connection between matrix functions and wavelets 

Definitions: A function, or a distribution, ip satisfying (2.3.7) is said to 

be refinable, the equation (2.3.7) is called the refinem.ent equation, or also, 
as noted above, the "scaling identity", and ip is called the scaling function. 
The coefficients a„ of (2.3.7) are called the masking coefficients. 

We will mainly concentrate on the case when the set {a„} is finite. But 
in general, a function ip G L'^ (M) is said to be refinable with scale number N 
if (p (x/N) is in the L^-closed linear span of the translates {(p> {x — fc)}^.^^ 
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L2 (R); see, e.g., [HSS96, SSZ99, StZh98, StZhOl]. 

Since there are refinement operations which are more general than scal- 
ing (see for example [DLLPOl]), there are variations of (2.3.7) which are 
correspondingly more general, with regard to both the refinement steps 
that are used and the dimension of the spaces. The term "scaling identity" 
is usually, but not always, reserved for (2.3.7), while more general refine- 
ments lead to "refinement equations". However, (2.3.7) often goes under 
both names. The vector versions of the identities get the prefix "multi-", 
for example multiscaling and multiwavelet. 

If mo satisfies a condition for obtaining orthogonal wavelets, 

J2 \mo{w)f = N, (2.3.42) 

together with the normalization 

mo (1) - V^, (2.3.43) 

then (2.3.7) has a solution (f in (R) which can be obtained by taking the 
inverse Fourier transform of the product expansion 

(Here and later we use the convention that if m (z) is a function of ^ e T, 
then m{t) =m (e"'*).) That (2.3.44) gives a solution ip of (2.3.7) follows 
from the relation 

2.3.3.1 Multiresolution wavelets 

We mentioned that there is a direct connection between mo = ^'^d 
the scaling function on M given in (1.1.34), (2.3.7), and (2.3.44). There is 
a similar correspondence between the high-pass filters m^ and the wavelet 
generators V'i S L'^ (K). In the biorthogonal case, there is a second system 
ifii ^ ipi and the two systems 

{Niipi {N^x - fc)} and {iV^^i/ (n^' x - fc')} , 

i, e {1, 2, . . . , AT _ 1} , j', fc, fc' e Z, (2.3.46) 

then form a dual wavelet basis, or dual wavelet frame for (M) in the sense 
of [Dau92], Chapter 5. We considered this biorthogonal case in more detail 
in § 2.3.1 above. Much more detail can be found in Chapter 6 of [BrJo02b]. 
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The idea of constructing maximally smooth wavelets when some side 
conditions are specified has been central to much of the activity in wavelet 
analysis and its applications since the mid-1980's. As a supplement to 
[Dau92], the survey article [Stra93] is enjoyable reading. The paper [LaHc96] 
treats the issue in a more specialized setting and is focussed on the moment 
method. Some of the early applications to data compression and image 
coding are done very nicely in [HSS + 95], [SHS + 99], and [HSW95]. An 
interesting, related but different, algebraic and geometric approach to the 
problem is offered in [PeWi99]. 

We now turn to an interesting variation of this setup, which includes 
higher dimensions, i.e., when the Hilbcrt space is (W^), d = 2,3, ... . 
Staying for the moment with d= 1, and N fixed, we will take the viewpoint 
of what is called resolutions, but here understood in a broad sense of closed 
subspaces: A closed linear subspace V C (M) is said to be an A''-resolution 
if it is invariant under the unitary operator 

U = Un: f^N-'^f(^^), (2.3.47) 

i.e., if U maps V into a proper subspace of itself. The subspace V is said to 
be translation invariant if 

/ G V / ( • - fc) G V for all fc e Z. (2.3.48) 

If there is a function (p such that V = V^ is the closed linear span of 

{ip{ - -k)\keZ}, (2.3.49) 

then clearly V is translation invariant. The translation- invariant resolution 
subspaces V are actively studied and reasonably well understood. If V is of 
the form in (2.3.49), then we say that it is singly generated, and that ip 
is a scaling function of scale N. 

2.3.3.2 Generalized multiresolutions [joint work with L. Baggett, 
K. Merrill, and J. Packer] 

The case when the resolution subspace V is not singly generated is also 
interesting, and these resolution subspaces are frequently called general- 
ized, multiresolution subspaces (GMRA). There is much current and very 
active research on them; see, for example, [BaLa99], [LPTOl], [BaMe99], 
[HLPS99], [HSSOl], [SSZ99], and [JorOla]. The case when V is not singly 
generated as a resolution subspace of scale N > 2, i.e., when V is not of 
the form (2.3.49), occurs in the study of wavelet sets. A wavelet set in M."^ is 
defined relative to an expansive dx d matrix N over Z. A subset E cM.'^is 
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said to be an N- wavelet set if there is a single wavelet function ip £ (W^) 
such that tp = xe- Specifically, the condition states that the family 

{ |det Nf^ ip (N^x - k) : j e Z, k e Z'' } (2.3.50) 

is an orthonormal basis for (M**) . This can be checked to be equivalent 
to the combined set of two tiling properties for E as a subset of Mf^: 

(a) the family of subsets { WE : j e Z} tiles R'^; 

(b) the translates {E + 27rk : k e Z'' } tile R''. 

We define tiling by the requirement that the sets in the family have over- 
lap at most of measure zero relative to Lebesgue measure on R"^. Similarly, 
the union 

K'^ = y -NiE = IJ £ + 27rk (2.3.51) 

jez kez<i 

is understood to be only up to measure zero. 

It is easy to see that compactly supported wavelets in (M.'^) are MRA 
wavelets, while most wavelets ip = (xe)^ from wavelet sets E are not. Thess 
wavelets are typically (but not always) frequency localized. 

The main difference between the GMRA (stands for generalized mul- 
tiresolution analysis) wavelets and the more traditional MRA ones may be 
understood in terms of multiplicity. Both come from a fixed resolution sub- 
space Vo C L"^ (M'') which is invariant under the translations { T„ : n e Z'* } 
where 

(T„/) (x) := f {x - n) for x G R"^ and n € Z'^. (2.3.52) 

Hence {IralvolngZ'^ unitary representation of Z'' on the Hilbert space 
Vo- As a result of Stone's theorem, we find that there are subsets 

El D E2 D ■ ■ ■ D Ej D ■ ■ ■ 

of T'' such that the spectral measure of the (restricted) representation has 

multiplicity > j on the subset Ej, j = 1,2, It can be checked that 

the projection- valued spectral measure is absolutely continuous. Moreover, 
there is an intertwining unitary operator 

J: Vo — > (^j) (2.3.53) 

such that 



PL2iE,)JTnf{z)=z"{Jf) (z) 



(2.3.54) 
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holds for all / <E Vq and z £ Ej. We may then consider the functions ipj G Vq 
(C (R'^)) defined by 

iPj := J-i(0,...,0, XB^,0,0,...). (2.3.55) 

j'th place 

It was proved by Baggett and Merrill [BaMe99] that {(pj : j > 1} generates 
a normalized tight frame for Vq: specifically, that 

|(T„<p, |/)|2 = ||/||i.(j,.) (2.3.56) 

holds for all / S Vq. 

Treating {ipi, (p2, ■ ■ ■) as a vector- valued function, denoted simply by (f, 
we see that there is a matrix function 



H : T"^ — > (complex square matrices) 



such that 



for transpose 

if (N*^f) = H (e'') ^ (t) , (2.3.57) 

where t={ti,...,td) G K'', and e^* := (e**i , e^*^ . . . , e^*-*) . 

But this method takes the Hilbert space (M'') as its starting point, 
and then proceeds to the construction of wavelet filters in the form (2.3.57). 
Our current joint work with Baggett, Merrill, and Packer reverses this. It 
begins with a matrix function H defined on T*^, and then offers subband con- 
ditions on the matrix function which allow the construction of a GMRA for 

(R'^) with generator ip = {cpi, ip2, ■ ■ ■) given by (2.3.57). So the Hilbert 
space (M.'^) shows up only at the end of the construction, in the conclu- 
sions of the theorems. 

2.3.4. Matrix completion 

In using the polyphase matrices, one may only have the first few rows, and 
be faced with the problem of completing to get the entire function A from 
a torus into the matrices of the desired size. The case when only the first 
row is given, say corresponding to a specified low-pass filter, is treated in 
[Br.Jo02b] and [BrJo02a], and we refer the reader to the references given 
there, especially [JiSh94], [RiShQl], [RiSh92], and [Vai93]. 

The wavelet transfer operator is used in a variety of wavelet applica- 
tions not covered here, or only touched upon tangentially: stability of refin- 
ablc functions, regularity, approximation order, unitary matrix extension 
principles, to mention only a few. The reader is referred to the following 
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papers for more details on these subjects: [DHRS03], [RoSh03], [RSTOl], 
[JJSOl], [RoShOO], [RiShOO], [JiSh99], [She98], [RoSh98], [LLS98], [RoSh97], 
[BJMP03], and [BJMP04]. 

The unitary extension principle (UEP) of [DHRS03] involves the inter- 
play between a finite set of filters (functions on M/Z), and a corresponding 
tight frame (alias Parseval frame) in (W^) . 

For the sake of illustration, let us take d = 1, and scaling number A'' = 
2, i.e., the case of dyadic framelets. Naturally, the notion of tight frame 
is weaker than that of an orthonormal basis (ONB), and it is shown in 
[DHRS03] that when a system of wavelet filters m,, i = 0, 1, . . . , r is given 
(mo must be low-pass), then the orthogonality condition on the rn,,;'s always 
gets us a framelet in (M), i.e., the functions corresponding to the high- 
pass filters nii, i — 1, . . . ,r, generate a tight frame for (M), also called a 
framelet. The correspondence rrii to ■^j is called the UEP in [DHRS03]. 

The orthogonality condition for rrii, i = 0, 1, . . . ,r, referred to in the 
UEP is simply this: Form an (r -|- l)-by-2 matrix-valued function F (x) 
by using m, (x), i = 0, 1, . . . ,r in the first column, and the translates of 
the mi's by a half period, i.e., {x + 1/2), i = 0, 1, . . . , r in the second. 
The condition on this matrix function F (x) is that the two columns are 
orthogonal and have unit norm in l"^ for all x. Note that we still get the 
unitary matrix functions acting on these systems, in the way we outlined 
above. But there is redundancy as the unitary matrices are (r + l)-by- 
{r + 1). The reader is refered to [DHRS03] for further details. 

We emphasize that several of these, and other related topics, invite the 
kind of probabilistic tools that we have stressed here. But a more systematic 
discussion is outside the scope of this brief set of notes. We only hope to 
offer a modest introduction to a variety of more specialized topics. 

Remark 2.3.4.1: The orthogonality condition for nij, i = 0, 1, . . . , r, may 
be stated in terms of the operators Si from equation (2.9), N = 2. For 
each i = 0, 1, . . . , r, define an operator on (R/Z) as in (2.9). Then the 
arguments from Section 2 show that the orthogonality condition for TOj, 
i = 0, 1, . . . , r, i.e., the UEP condition, is equivalent to the operator identity 
(2.8) where the summation now runs from to r. Operator systems 5, 
satisfying (2.8) are called row-isometries. 

Remeirk 2.3.4.2: There are two properties of the low-pass filter mo which 

wc have glossed over. First, toq must be such that the corresponding scaling 
fmiction ip is in (K) . Without an added condition on mo , might only be 
a distribution. Secondly, when the dyadic scaling in (M) is restricted to 
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the resolution siibspacc Vq {(p), the corresponding unitary part must be zero. 
These two issues are addressed in [BJMP03], [BJMP04], and [DHRS03]. 

2.3.5. Connections between matrix functions and signal processing 

Since our joint work with Baggett, Merrill, and Packer on the GMRA 
wavelets is still in progress, we restrict the discussion of matrix functions 
here to the MRA case. 

The two groups of matrix functions C (T, Uat (C)) and C (T, GLjv (C)), 
i.e., the continuous functions from the torus into the respective groups, 
enter wavelet analysis via the associated wavelet filters {mi)f^^ . 

In [BrJo02b] (see also § 1.1.3 above), we give the details of the multiple 
correspondence between: 

(i) matrix functions. A: T ^ GLn (C), 

(ii) high- and low-pass wavelet filters mi, rhi', i,i' = 0,1, ... ,N — 1, and 

(iii) wavelet generators tpi, tpi', i,i' = 1, . . . ,N — 1, together with scaling 
functions (p, (p. 

In particular. 



The dependence of the 1? (IR)-functions in (iii) on the group elements 
A from (i) gives rise to homotopy properties. The standard orthogo- 
nal wavelets represent the special case when = m,, or equivalently, 
A {z) = ((A (-z))*) , 2 G T. Hence, the matrix functions are unitary in 
this case. 

The scaling/ wavelet functions V'l, • ■ • , V'iv-i with support on a fixed 
compact interval, say [0, kN -h 1], A; = 0, 1, . . . , can be parameterized with 
a finite number of parameters since the unitary- valued function z ^ A{z) 
in (2.3.58) then is a polynomial in z of degree at most k{N —1). It is 

well-known folklore from computer-generated pictures that the shape of 
the scaling/ wavelet functions depends continuously on these parameters; 
see Figures 1.1-1.7 in [BrJo02b] and [TreOl]. 

The scaling function ip G L'^ (K) of (2.3.7) is illustrated there, in the 
case TV = 2, and for orthogonal Z-translates, i.e., the case (2.3.42). These 
pictures illustrate the dependence of ip on the masking coefficients (o„) in 




(2.3.58) 



(2.3.59) 
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the case of [TreOl]: 

flo = ivo -rii-V2 + V3+ '?4)/4, 
ai = (770 + »7i - »72 + ??3 - ??4)/4, 
(12 = {vo -m- V4:)/2, 
as = {rjo -V3+ Vi)/^, 
ai = ivo + Vi + m + V3 + »?4)/4, 
as = {vo -111+112 + 113- '?4)/4, 

where 

Vo = l/\/2, 

r?i = (cos26i + cos2/9)/\/2, r?2 = (sin26i + sin2p)/\/2, (2.3.61) 

% = cos(26' - 2p) /V2, r]4 = sm{2e - 2p) /\/2. 

These formulas arise from an independent pair of rotations by angles 9 and 
p of two "spin vectors", i.e., by taking the matrix function A in (2.3.58) 
unitary, T 9 2; — > Ae^p (2) € U2 (C), and setting 

A{z) = V{Qj + zQg){Q^ + zQp) = VU6{z)Up{z) {'ZJim) 

with 

Q _ f cos^ 6 cos 6 sin 6 
\ cos 9 sin 9 sin^ 9 

lfflQ\ f cos 29 sin29 \\ /^o^.n 
= 2(loij + (sin2^-cos2^jj' (^'^-^^^ 

and the orthogonal complement to the one- dimensional projection Qg, 

Qi = Qe+(./2y (2-3.65) 

With the coefficients og, oi, 02, 03, 04, as given by (2.3.60), the algorith- 
mic approach to graphing the solution to the scaling identity (2.3.7) is as 
follows (sec [JorOlb], [TrcOl] for details): the relation (2.3.7) for W = 2 is in- 
terpreted as giving the values of the left-hand (p by an operation performed 
on those of the ip on the right, and a binary digit inversion transforms this 
into the form 

K+iU+^)=^fkix), (2.3.66) 
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where A is the 2x3 matrix Aj j = \/2a4+i-2i constructed from the coeffi- 
cients in (2.3.7), and fj and f are the vector functions 




(2.3.67) 



The signal processing aspect can be understood from the description 
of subband fihers in the analusis and synthesis of time signals, or more 
general signals for images. In either case, we have two subband systems 
m = (mo, nil, • ■ • ) and rh = (mo, m-i, . . . ) where the functions 



arc the generating functions defined from the filter coefficients ( a,'' ) and 



Appendix A: Topics for further research 

Originally we had anticipated adding two more chapters to these tutorials, 
but time and space prevented this. Instead we include the table of contents 

for this additional material. The details for the remaining chapters will be 
published elsewhere. But as the items in the list of contents suggest, there 
are still many exciting open problems in the subject that the reader may 
wish to pursue on his/her own. We feel that the following list of topics 

offers at least an outline of several directions that the rc^adcr, could take in 
his/her own study and research on wavelet-related mathematics. 



m. 



.,-(^)=^a(^)^" and m,(z)=^aW^^ 



n 



n 
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3. Connection between the discrete signals and the wavelets 

3.1. Wavelet geometry in L^{R'^) 

3.2. Intertwining operators between sequence spaces P and 



3.3. Infinite products of matrix functions 

3.3.1. Implications for L'^{R'^) 

3.3.2. Wavelets in other Hilbert spaces of fractal measures 

3.4. Dependence of the wavelet functions on the matrix 
functions which define the wavelet filters 

3.4.1. Cycles 

3.4.2. The Ruelle-Lawton wavelet transfer operator 
4. Other topics in wavelets theory 

4.1. Invariants 

4.1.1. Invariants for wavelets: Global theory 

4.1.2. Invariants for wavelet filters: Local theory 

4.2. Function classes 

4.2.1. Function classes for wavelets 

4.2.2. Function classes for filters 

4.3. Wavelet sets 

4.4. Spectral pairs 

Appendix B: Duality principles in analysis 

Several versions of spectral duality are presented. On the two sides we 
present (1) a basis condition, with the basis functions indexed by a fre- 
quency variable, and giving an orthonornial basis; and (2) a geometric 
notion which takes the form of a tiling, or a Iterated Function System 
(IFS). Our initial motivation derives from the Fuglcdc conjecture, sec 
[Fug74, Jor82, JoPe92]: For a subset D of M" of finite positive measure, 
the Hilbert space {D) admits an orthonormal basis of complex exponen- 
tials, i.e., D admits a Fourier basis with some frequencies L from M", if 
and only if D tiles R" (in the measurable category) where the tiling uses 
only a set T of vectors in M". If some D has a Fourier basis indexed by a 
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set L, we say that {D, L) is a spectral pair. We recall from [JoPe99] that 
if D is an n-cube, then the sets L in (1) are precisely the sets T in (2). 
This begins with work of Jorgensen and Steen Pedersen [JoPe99] where the 
admissible sets L = T are characterized. Later it was shown, [IoPe98] and 
[LRWOO] that the identity T = L holds for all n. The proofs are based on 
general Fourier duality, but they do not reveal the nature of this common 
set L = T. A complete list is known only for n = 1, 2, and 3, see [JoPe99]. 

We then turn to the scaling IFS's built from the n-cube with a given 
expansive integral matrix A. Each A gives rise to a fractal in the small, and a 
dual discrete iteration in the large. In a different paper [JoPe98], Jorgensen 
and Pedersen characterize those IPS fractal limits which admit Fourier 
duality. The surprise is that there is a rich class of fractals that do have 
Fourier duality, but the middle third Cantor set does not. We say that an 
afRne IFS, built on afHne maps in M" defined by a given expansive integral 
matrix A and a finite set of translation vectors, admits Fourier duality if the 
set of points L, arising from the iteration of the A-affine maps in the large, 
forms an orthonormal Fourier basis (ONB) for the corresponding fractal /i 
in the small, i.e., for the iteration limit built using the inverse contractive 
maps, i.e., iterations of the dual affine system on the inverse matrix A~^. By 
"fractal in the small" , we mean the Hutchinson measure fj, and its compact 
support, see [HutSl]. (The best known example of this is the middle-third 
Cantor set, and the measure fj, whose distribution function is corresponding 
Devil's staircase.) 

In other words, the condition is that the complex exponentials indexed 

by L form an ONB for L'^ifJ,)- Such duality systems are indexed by complex 
Hadamard matrices H, sec [JoPc99] and [JoPe98]; and the duality issue is 
connected to the spectral theory of an associated Ruelle transfer operator, 
see [BrJo02b]. These matrices H are the same Hadamard matrices which in- 
dex a certain family of quasiperiodic spectral pairs (D, L) studied in [Jor82] 
and [JoPe92]. They also are used in a recent construction of Terence Tao 
[Tao04] of a Euclidean spectral pair {D, L) in for which D does not a 
tile with any set of translation vectors T in R^; see also [IKT03]. 

We finally report on joint research with Dorin Dutkay [DuJoOS], 
[DuJo04a], [DuJo04b], [DuJo04c] where we show that ah the affine IFS's, 
and more general limit systems from dynamics and probability theory, ad- 
mit wavelet constructions, i.e., admit orthonormal bases of wavelet func- 
tions in Hilbert spaces which are constructed directly from the geometric 
data. A substantial part of the picture involves the construction of limit 
sets and limit measures, a part of geometric measure theory. 
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